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Final  Progress  Report  -  DAAD19-01-1-0631 

Purpose  and  Goals  of  Effort 

The  overarching  goals  of  this  effort  has  been  the  continued  development  of  high-order 
accurate  computational  methods  for  wave  problems  in  general  and  Maxwell’s  equations 
of  electromagnetics  in  particular.  The  applications  are  general  but  particular  emphasis 
has  been  on  optical  applications. 

Rather  than  focusing  on  one  particular  approach,  we  have  pursued  such  developments 
along  three  different  directions 

•  Time-Domain  High-order  Embedding  Methods  Focuses  on  the  development,  analy¬ 
sis,  and  implementation  of  a  4th  order  finite  difference  embedded  interface  scheme 
for  the  time-domain  solution  of  Maxwell’s  equations.  This  enables  the  solution  of 
geometrically  complex  problems,  including  real  materials,  using  a  Cartesian  grid 
without  sacrificing  accuracy,  stability  and  robustness.  Part  of  this  also  requires  at¬ 
tention  to  related  issues  such  a  efficient  and  robust  absorbing  boundary  conditions. 

•  Boundary  Variation  Methods  for  Multilayered  Optics  Focuses  on  boundary  varia¬ 
tion  techniques  for  the  solution  of  Maxwell’s  equations  in  the  frequency  domain. 
Applications  are  taken  from  photonics  and  integrated  optics,  including  waveguide 
couplers,  multilayer  devises,  and  transmission  and  reflection  diffractive  optics. 

•  DG-FEM  Methods  for  the  Harmonic  Maxwell  Problem  Focuses  on  the  development 
and  analysis  for  Discontinuous  Galerkin  Methods  for  the  harmonic  frequency  do¬ 
main  Maxwell  problem  as  well  as  the  wave  Helmholtz  problem.  These  are  generic 
methods,  the  details  of  which  are  still  poorly  understood  and  we  have  developed  the 
analysis  and  computational  tools  in  parallel.  Subtle  issues  such  as  spuriuos  modes 
and  efficient  preconditioning  techniques  for  the  resulting  indefinite  linear  systems. 

Overview  of  Main  Results 

In  the  following  we  shall  discuss  these  efforts  in  some  more  detail,  highligthing  the 
progress  made  during  this  effort  as  reported  in  the  publications. 

Embedded  Interface  Methods 

Within  the  scope  of  this  research,  we  have  achieved  the  following  most  notable  tasks 

•  Following  previous  pioneering  work  by  us,  we  have  developed  4th  order  accurate 
finite  difference  schemes  which  are  able  to  solve  wave  problems  stated  in  complex 
geometries.  This  development  overcomes  classic  problems  of  staircasing  and  the 
effect  internal  boundaries  in  a  unified  way.  As  the  scheme  is  modified  locally  only, 
it  maintains  the  simplicity  and  computational  efficiency  of  the  original  formulation 
as  most  of  the  additional  work,  i.e.,  computing  the  local  schemes,  is  done  in  a 
preprocessing  stage. 
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The  extension  of  these  ideas  to  4th  order  embedding  methods  is  far  from  trivial. 
However,  we  have  succeeded  in  formulating  and  implementing  such  a  scheme,  using 
a  standard  staggered  grid  everywhere  in  space  and  locally  modifying  the  scheme  by 
use  of  one  sided  stencils  and  extrapolations. 

Computational  experimentation  and  preliminary  analysis  confirms  the  stability  of 
the  embedding  scheme  for  situations  where  the  material  interfaces  are  straight.  The 
scheme  is  currently  being  implemented  and  tested  in  a  two-dimensional  version, 
with  preliminary  results  being  very  encouraging. 

The  scheme  is  currently  being  implemented  in  a  production  code  for  use  on  DoD- 
related  applications. 

•  We  have  investigated  some  long  time  stability  problems  observed  in  PML  methods 
used  for  absorbing  waves  in  the  time-domain  schemes.  Extensive  computations 
showed  the  potential  for  weak  instability  and  the  source  of  this  was  found  to  be 
in  the  degenerate  low-order  terms.  This  understanding  enabled  stabilization  by 
breaking  the  degeneracy. 

•  We  have  proposed  an  entirely  new  family  of  absorbing  layer  methods  of  a  nonlinear 
nature  -  known  as  nonlinear  PML  methods.  These  are  perfectly  matched  layers, 
yet  require  no  additional  variables,  thus  reducing  the  memory  usage  signifcantly. 
The  analysis  of  this  scheme  remains  incomplete  due  to  its  nonlinear  nature  but 
computational  results  very  encouarging. 

Boundary  Variation  Methods 

Within  the  scope  of  this  research,  we  have  achieved  the  following  most  notable  tasks 

•  A  first  goal  was  to  finalize  work  on  the  development  of  three-dimensional  boundary 
variation  methods  for  single  interface  scattering,  in  particular  with  applications  to 
waveguide  couplers. 

This  has  been  accomplished  as  planned  and  several  large  scale  applications  has  been 
completed. 

•  As  a  second  goal,  we  have  developed  an  entirely  new  class  of  boundary  variation 
methods  for  multilayered  transmission  optics. 

A  flexible  and  fast  approach  for  the  two-dimensional  harmonic  problem  has  been 
developed.  The  algorithm  relies  on  recursive  use  of  a  highly  optimized  one-layer 
scheme  and  allows  the  modeling  of  an  arbitrary  number  of  dielectric  layers,  each 
with  an  interface  of  arbitrary  shape.  Important  special  cases  include  the  ability  to 
model  transmission  optics  of  finite  extend  and  problems  illuminated  by  beams. 

A  key  development  is  the  use  of  physical  insight  to  enable  a  fast  and  cost  effective 
of  the  recursive  approach.  With  out  relying  on  the  physical  properties  of  the  scat¬ 
tering  and  diffraction  processes,  the  recursive  process  would  grow  exponentially  in 
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complexity.  With  this  new  formulation,  the  growth  is  bounded  from  above  and  all 
central  components  are  precomputed. 

The  developed  software  has  been  validated  carefully  by  comparison  with  selected 
exact  solutions  as  well  as  by  careful  cross  validation  of  results  obtained  using  the 
high-order  time-domain  codes.  Extensions  to  the  three-dimensional  vectorial  case 
is  expected  to  be  relatively  straightforward,  based  on  the  general  algorithmic  devel¬ 
opments. 

•  The  methods  for  fast  and  accurate  modeling  of  multilayered  diffractive  optics  have 
been  validated  and  used  extensive  by  several  research  groups,  e.g.,  at  Bell  Labs  for 
studies  of  the  impact  of  roughness  in  Bragg  mirrors  in  lasers  and  at  the  National 
Lab  of  Denmark  to  model  and  design  highly  sensitive  sensors  for  measuring  bacterial 
content  in  liquids.  These  collaborations  have  both  led  to  further  developments  of 
the  algorithm  as  well  as  facilitated  an  increased  understanding  and  analysis  of  the 
particular  applications  at  hand.  The  results  of  both  efforts  are  being  prepared  for 
publication. 

•  In  a  closely  related  effort  we  have  demonstrated  the  use  of  the  boundary  variation 
method  as  a  tool  to  analyze  errors  coming  with  poor  geometry-descriptions  in  stan¬ 
dard  finite-difference  schemes  used  for  solving  Maxwell’s  equations.  The  boundary 
variation  method  allows  for  solving  the  diffraction  problem,  appearing  by  represent¬ 
ing  a  straight  interface  by  a  stair-cased  approximation,  exactly,  hence  providing  a 
mechanism  to  separate  discretization  errors  from  errors  from  geometry  descriptions. 
The  results  show  that  the  physics  of  the  problems  changes  dramatically  and  one  can 
only  expect  first  order  convergence  in  amplitude  and  phase  computed  fields.  This 
is  indeed  a  strong  argument  for  why  embedded  interface  methods  are  needed. 

DG-FEM  Methods  for  the  Harmonic  Maxwell  Problem 

Within  the  scope  of  this  research,  we  have  achieved  the  following  most  notable  tasks 

•  we  We  have  initiated  the  investigation  of  the  ability  to  accurately  and  efficiently  ap¬ 
proximate  eigenvalue  problems  originating  in  Maxwell’s  equations,  i.e.,  when  com¬ 
puting  eigenmodes  of  resonators.  This  is  a  notoriously  hard  problem  due  to  the 
nature  of  the  curl-curl  form  of  Maxwell’s  equations  and  the  associated  boundary 
conditions. 

The  emphasis  has  been  on  the  possibility  of  accurately  and  reliably  solve  the  eigen¬ 
value  problem  for  nonsmooth  geometries  where  the  eigensolutions  can  exhibit  strong 
singularities. 

The  computional  approach  is  that  of  high-order  discontinuous  Galerkin  methods  on 
fully  unstructured  grid.  Such  methods  has  not  been  explored  in  detail  for  Maxwell 
eigenvalue  problems  previously  and  our  study  analysis  has  revealed  some  unexpected 
properties  of  the  discretization,  the  discrete  null  space  of  the  operators  and  ways  to 
stabilize  it. 
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The  study  has  been  focused  on  the  possibility  of  exciting  spurious  modes,  as  well- 
known  problem  in  classical  nodal  finite  element  methods,  and  general  behavior.  No 
significant  problems  of  spurious  modes  have  been  found,  confirming  that  this  formu¬ 
lation  is  well  suited  for  dealing  with  geometrically  complex  nonsmooth  problems. 

•  In  the  frequency  domain  we  have  continued  the  development  of  high-order  accu¬ 
rate  discontinuous  Galerkin  methods  for  Maxwell’s  equations.  The  schemes  are 
completely  general  and  highly  flexible,  using  arbitrary  order  and  fully  unstructured 
grids. 

Solving  the  resulting  discrete  problems  is  challending  for  several  reasons.  The  size 
can  very  considerable  and  we  must  seek  techniques  that  map  to  parallel  computing 
platforms.  Furthermore,  the  condition  number  of  the  matrix  scales  quite  poorly  with 
both  the  inverse  element  size  and  the  order  of  approximation.  Finally,  the  problem 
is  highly  indefinite  which  rules  out  most  standard  preconditioning  techniques. 

As  a  prototype  problem  we  are  considering  the  wave  Helmholtz  problem  for  which  we 
have  developed  a  two-level  additive  overlapping  Schwarz  preconditioner  and  tested 
it  extensively.  In  the  low  frequency  case,  where  the  matrix  remains  definite,  this  pre¬ 
conditioner  complete  removes  both  grid  and  order  dependence,  resulting  in  constant 
iteration  count. 

For  the  indefinite  case,  the  coarse  grid  solver  becomes  essential  in  reducing  the  de¬ 
pendence  of  the  iteration  count  on  the  wavenumber.  However,  it  does  not  eliminate 
it  completely  as  eigenvalues  cross  the  axis  when  being  shifted  from  coarse  to  fine 
grid  and  back.  The  number  of  these  grows  like  the  squareroot  of  the  wavenumber 
and  the  iteration  count  scales  like  this. 

However,  this  weak  dependence  is  better  that  most  other  schemes  and  the  over¬ 
lapping  Schwarz  scheme  is  superior  when  also  taking  into  account  its  successful 
treatment  of  by  grid  and  order  dependence  as  well  as  it  inherent  parallel  nature. 

Currently  we  are  extenting  these  ideas  to  the  full  curl-curl  form  of  Maxwell’s  equa¬ 
tions  as  well  as  considering  the  optimization  of  fluxes  and  use  of  multiplication 
Schwarz  methods. 
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